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Abstract 

We investigate the modification of the pion-cloud in the nuclear medium and 
its effect on the nuclear Drell-Yan process. The pion's in-medium self-energy 
is calculated in a self-consistent delta-hole model, with particle-hole contribu- 
tion also included. Both the imaginary and real part of the pion's and delta's 
self-energy are taken into account and related through a dispersion relation 
assuring causality. The resulting in-medium pion light-cone momentum dis- 
tribution shows only a slight enhancement compared to the one of the free 
nucleon. As a consequence the ratio of the cross-section for Drell-Yan scat- 
tering on nuclear matter and nucleonic target is close to unity in agreement 
with experiment. 
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I. INTRODUCTION 



Recently there has been renewed interest in the role of the pion propagator in the nuclear 
medium for several reasons: First the ratio of spin-longitudinal and spin-transverse response 
functions below the quasi-elastic peak in a naive model is predicted to be much larger than 
unity, while experiment finds a ratio close to unity in (p, n) polarization transfer [JTJ] . Second 
the observed ratio of cross sections for Drell-Yan scattering on nuclear and nucleonic targets 
is consistent with no excess pions present in the nucleus, i.e. no enhancement of the sea- 
quarks p. 

In the past the effect of medium modifications of the pion propagator has been in- 
vestigated in detail by Ericson and Thomas || (in connection with the EMC effect) and 
Bickerstaff et al. [Q . The apparent absence of the pion-cloud enhancement in DY scattering 
was explained by Brown et al. || in terms of partial restoration of chiral symmetry and the 
associated decrease in masses of the nucleons and vector mesons in the nuclear medium. It 
was also pointed out that, in a definite model considered, the correct normalization of 
the physical nucleon state, consisting of the bare nucleon and pion-cloud term, considerably 
reduces the effect of the pion-cloud enhancement. 

The sensitivity of the pion-cloud to the pion-nucleon-nucleon (ttNN) vertex cut-off A was 
emphasized by Thomas |7]]. Not only is the sea-quark content of the nucleon very sensitive 
to the value of the cut-off, the enhancement of the pion light-cone-momentum distribution in 
the nuclear medium also follows that pattern, even if the value of the cut-off is kept the same 
in the medium as in free space |§. Another uncertainty concerns the values of the Migdal 
g' parameters (s'tvtv 9na^ 9aa)j describing the effects of short-range repulsion. Taking large 
enough values of these parameters it is possible to reduce the medium enhancement of the 
pion distribution. However, it has not been reliably established that using such large values 
of the g' parameters is physically justified. 

Our approach in reexamining the pion-cloud enhancement is to perform a more ambitious 
computation of the pion light-cone-momentum distribution in the nuclear medium, using the 
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energy and momentum dependent pion self-energy and propagator. The latter are calculated 
in a recently developed self-consistent delta-hole model, including particle-hole states, and 
allowing for a width of the delta |J. This model takes into account both the real and the 
imaginary parts of the pion's and delta's self-energy, in a way that assures causality and 
absence of unphysical poles. The Schwinger-Dyson equations (without vertex corrections) 
for the delta and the pion are solved self-consistently, while the nucleon is treated in a 
mean-field approximation. 

The framework of the present approach (and that of Ref. |J) is the effective field theory of 
hadrons, which is used consistently to calculate physical processes. The basic difference with 
previous approaches is the self-consistent treatment of the delta-hole contribution to the pion 
self-energy, including both its real and imaginary part. As already noted in Ref. || much 
of the pion-distribution enhancement came from the delta-hole term through its coupling to 
particle- hole term by g^A- However, that result was based on a crude approximation for the 
delta- hole contribution, including only its real part (violating causality), with a momentum 
dependence coming only from the p-wave nature of the coupling and the form factor. 

For the practical calculation of the pion light-cone momentum distribution in isospin 
symmetric nuclear matter (at temperature T = 0) we compare two approaches. First in 
section II we present a calculation of the pion light-cone-momentum distribution, based 
on explicit summation of the relevant diagrams involving the dressed pion propagator, and 
using an explicit integration over the initial and final nucleon momenta. 

Second, in section III we consider the computation of the same quantity, but based on 
the nuclear-matter response function. While somewhat simpler, this approach gives only 
an approximate (however, for typical nuclear densities very close) expression for the pion 
distribution. This happens since the calculation of the imaginary part of the pion self-energy 



involves the full phase-space factors for the two nucleons ||10|| , while such momentum de- 
pendent factor is not present for the incoming nucleon in the calculation involving explicit 
summation of relevant diagrams. However, the response-function approach takes into ac- 
count also the process when a nucleon becomes a delta plus a pion, if the pion self-energy 



includes the delta-hole contribution. We checked numerically that the two approaches give 
practically the same result. The ratio for Drell-Yan scattering on an isospin-symmetric nu- 
clear target and the deuteron is presented in section IV. Numerical results, a discussion and 
comparison with other computations and experimental results is the subject of section V. 



II. PION DISTRIBUTION IN THE NUCLEAR MEDIUM 



It was noticed long time ago that the pion cloud gives a scaling contribution to the 



deep inelastic lepton scattering on the free nucleon ||11|| . The nucleon wave function can be 
schematically expressed as 

| JV > ph ys= y/Z\N > bare +a|iVvr > +/3| Avr > + • • • (1) 

In this approach the light-cone momentum distribution of a quark with flavor / in a proton 
can be written as (B = N, A) 



q f (x) = Zq fMTe (x) + 



B.i 



1 -f B > /N (yhf,iare(x/y) + f dj -r ,N {v)^ (x/y) 

x y Jx y 



(2) 



Here Cj (i labels the charge states) are the appropriate isospin Clebsch-Gordan coefficients. 
Note that following |T2|] we use the Z factor, the bare nucleon probability, as a renormaliza- 



tion of the bare nucleon only, not as an overall normalization of the whole right-hand side. 
In Refs. |5]J| the Sullivan contribution was also multiplied with the wave function renor- 
malization factor. To see the relation we note that in quantum field theory one can write 
\N > p hys= \ // ~Z'{\N >barc +fl , °Arvl^ 7r >)> where g° denotes the bare coupling. The latter is 
related to the physical one (in lowest order) by g = y/Z'g , i.e. Z' = (1 + / dyf n / N (y))-\ 



The prescription used in the present paper (see also Ref. fl3f ) is consistent with the stan- 
dard nuclear physics definition of the tcNN coupling constant derived from NN interaction 
at large distances. Of course the basic assumption is that we can restrict ourselves to only 
one pion in the air. For the form-factors we used in present calculation (whose conclusions 
are not affected by some variation of the cut-off) the value of Z 0.6 — 0.7 suggests that 
this is a good approximation. 



The pion light-cone momentum distribution, f w ' N (y), may contain either a nucleon or a 
delta final state i.e. f n ^ N (y) = f nN ^ N (y) + f* A ^ N {y)- Let us first consider the nuclear final 
state. For the free nucleon the well-known result for the pion distribution is 



r w [y) = 3 Am y r dt 'Wf , (3) 



167T 2 JM 2 y 2 /{l-y) (t + W 2 ) 2 

Here y = k °~j^ 3 gives the pion light-cone momentum fraction, with M being the physical 
mass of the nucleon (as an arbitrary but convenient scale), F n NN(t) is the nNN vertex 
form-factor, while g n NN is the n°NN coupling. The free-pion propagator, _D°, appears in 
the above expression in the form (t + m^) _1 . 

In the literature various prescriptions for the irNN form-factor have been proposed. 
Instead of the covariant formalism, with the pion taken off-shell and the final nucleon on- 
shell as used above, in Ref. |L2| and also in work of the Jiilich group, Ref. ||, old-fashioned 



perturbation theory in the infinite momentum frame was used. In that formulation the pion 
and final nucleon are taken to be on the mass shell (however, energy is not conserved in the 
interaction vertex) and the form-factor is conveniently taken a function of the center-of-mass 
energy of the intermediate ttN system, s = (p' + k) 2 . In this case the probability to find a 
baryon in nucleon with momentum fraction y is equal to the probabilty of finding a meson in 
the nucleon with momentum fraction 1 — y, which is not the case in the covariant approach 
with the form-factor depending only on t = —k 2 . 

The impracticality of making in-medium (especially self-consistent) calculations in the 
infinite-momentum frame makes us adopt the covariant approach. Since the covariant and 
the old-fashioned perturbation approach are equivalent (apart from the effect of form fac- 
tors), we assume the conditions which should be satisfied by the first and second moments 
of the two distributions for flavor-charge and momentum conservation can be assured (at 
least approximately), also in the former scheme by a suitable p' 2 dependence in the form- 
factor. In that case (see section IV) it is enough to consider only the change of f*{y) due 
to medium effects, which is the main subject of the present work. Furthermore, the pion 
distributions obtained in the two schemes are very similar if in the covariant approach a 



monopole form-factor is used as in our calculation. 



It was confirmed |TJ| that dressing the pion in vacuum gives practically negligible contri- 
bution. That is not the case, however, in nuclear medium where delta-hole and particle-hole 
states play an important role. Consequently, we want to base our calculation on the dressed- 
pion self-energy and propagator, depending on fc and k = \k\. The types of diagrams which 
contribute are shown in Fig. 1. They can be summed in RPA in two steps: 
(i) first, the pion rescattering (without the g' coupling to external nucleon) leads to the 
dressed pion propagator D n (k ,k), D^ 1 = (D®)^ 1 — k 2 U, where the pion self-energy k 2 U 



contains both the delta-hole and particle- hole contributions ||15|| : 



1 ~~ (,9na) N- NN U na 

Here ITvAr and IT at a are the iterated particle-hole and delta-hole self-energies divided by the 
common factor k 2 and the arguments ko, k are left out for brevity: 

rijvjv — n^jy^o, k) 

1 ~~ 9nn^°nn(^o^ k) 
t-t _ n^ A (fc , k) . . 

with H% N and H% A being the one-loop particle-hole and delta-hole self-energies devided by 
k 2 . 

(ii) second, the inclusion of diagrams in which the pion first couples with g' has the effect 
of replacing D.,, by 

n (v v\ n (u u\ (, , X 1 (k ,k)+X 2 (k ,k) \ 

D irN {k ,k) = D ir {ko,k)[l + -— — (7) 

\ X 3 (k ,k) J 

with 

Xx(ko, k) = g' NN IL NN (1 + g' NA IL NA ) , (8) 
X 2 (k , k) = g' NA n N A (1 + Qna^-nn) , (9) 

x 3 (k , k) = [i - {g' NA f n^n^] • (10) 

As a result in nuclear matter D® in Eq. (^) is replaced by and in addition the integration 
limits are adjusted. The nuclear matter ground state is approximated by a mean field, i.e. 
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the nucleon that emits the pion has an initial momentum p < pp and to satify the Pauli 
principle a final momentum p' > pp. The effect of binding is included in an effective mass 
and a shift of energy approximation, i.e. E = ^Jp 2 + M 2 + Cq. The values of the parameters 
in the above expression are determined from nuclear matter models and the requirement for 
leading to a reasonable value of the Fermi energy. 

As a result the pion light-cone distribution (per nucleon) in isospin symmetric nuclear 
matter can be expressed as 

x rd9k 2 F^ NN {k)^\D^ N \ 2 . (11) 
Jo z' 

In the above p'± max = ^jlz'M^Ep — (z' 2 M 2 + M 2 ), 9 is the angle between p± and p±, k = 
p — p' (the three-momentum of the pion) , and 

z ' = w, i~ My +P3 + +vl+ p2± ) = wX p ' z + + p ' 2 ) ' (12) 

ko = \[M* +Pl+Pl - ^M 2 +p>l+p>i, (13) 

We note that f n ^ A (y) has the correct non-relativistic low-density limit, \im PF ^of n ^ A (y) = 
f n ^ N (y), which is not the case in Ref. §1. 

We also consider diagrams with a delta replacing the nucleon in the final state, in which 
case the relativistic analog of expression ([|) reads 

where M± = Ma ± M. The above expression takes into account the width of the delta 
through its (vacuum) spectral-function p&. The energy of the delta is determined by the 
light-cone-momentum of the pion as p$ = M(l — y) — p 3 and t = p\ — M 2 y 2 — 2Myp 3 . 
We used the nonrelativistic limit of this expression, which amounts to replacing the term 
(M+ + t) 2 (M^ + t) by the first term of low-momentum expansion: lQM 3 Mj\k 2 , where k is 
the magnitude of the pion's three-momentum. In isospin-symmetric nuclear medium the 
pion light-cone distribution (per nucleon) takes a form analogous to fllll) : 

7 



where 



with 



x d6k 2 Fl N ^k)\p A { P ^ P >)\D nA \ 
jo p Q 



P±dp± / dp 3 / p^dpj 



D wA (k , k) 



D n (k , k) 



1 ~ (fl^A) 2n ^njvA ~[ 



^2X Ai (k ,k), 



(15) 



(16) 



X A1 (/c , k) = g' NA H-NN (1 + 9'na^na) 
X A2 (k , k) = g' AA Il NA (1 + g' NA Il NN ) 
X A3 (k , k) = 1 - (^ A ) 2 njvjvIIivA 



(17) 
(18) 
(19) 



In Eq. (f5j) p' = ^ + p 3 2 , p' = ^jMl + p 2 + c - My + p 3 - p' 3 and k = y Ml + p\ + p\ + 
c o ~~ Po- Obviously the total pion light-cone momentum distribution in the medium is the 
sum of the contributions with a nucleon and with a A in the final state, 



r A {y) = r N A {y) + fl A {y). 



(20) 



III. PION DISTRIBUTION AND THE NUCLEAR RESPONSE FUNCTION 

The spin-isospin nuclear response function essentially involves the imaginary part of the 
(iterated) pion self-energy and since diagrams of the same type appear in deep inelastic 
scattering off the pion emitted by the nucleon, the pion distribution in the medium can be 
expressed in terms of the response function [f|]. 

We want to clarify the relation of these two quantities for the relativistic approach. 
Even though we are using nonrelativistic nucleon and delta propagators (as well as nNN 
and 7riVA vertices), our use of relativistic kinematics necessitates the discussion, especially 
for establishing the correct low-density limit, as well as for getting the factors of nucleon's 
physical and effective mass in the expression correctly (the latter does not appear in previous 
treatments H,[l6[|, while in Ref. || a peculiar rescaling is introduced). 
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As a first step we consider the relation of the imaginary part of the particle-hole pion self- 
energy and the relevant part of the deep inelastic diagram shown in Fig. 2. Iterating the pion 
self-energy (with free-pion propagators in between) to get the full response function will not 
change the established relation (in the kinematic region of interest), since the off-shellness of 
the emitted pion means that cutting the free-pion propagator gives zero contribution. The 
square of the absolute value of the diagram shown in Fig. 2b gives: 

L = ~2gl NN I d3p /*T = TV[ 75 (A + + M*)]&(p F - 0)6(^1 - p P ) 

J (2vr) 6 2 J Ml + p' 2 



d 3 p d 3 k 
]M* + ( 

where the isospin degeneracy factor of 2 is included, p* denotes a p whose zeroth component 



-4£nn i / Pdk 5 (fcg - k 2 )Q(PF - \P\MW- k\ -Pf), (21) 

J (27r) 6 2yM 2 + (p — k) 2 



has a mean-field shift, 6 is the step function and k = y M 2 + p 2 — y M 2 + (jp — k) 2 . 

On the other hand, the contribution of the mean-field nucleons to the imaginary part of 



the pion self-energy in isospin symmetric nuclear matter can be written as fll7| , [14j 

Indigo, k) = -2 x 2gl NN n 2 J -^-Tr[ l5 (fa + M*) 75 ( A~ A + Af*)]<J((Po - c ) 2 - p 2 - M 2 
x5((p - c - k ) 2 - (p- k) 2 - Ml)Q{p F - \p\)Q{\p- k\-p F ) 

2 _ f (Pp 2 £*2\r//„ jl N2 /-» A2 7i^-2\ 



= -4g 2 NN 7r / " (fcg - k 2 )5((p ~ k ) 2 - (p- k) - M, 

J (2^2^ M 2 +p 2 

xQ{p F -\p\)e{\p-k\- PF ). (22) 

Integrating the above expression over fc an d k we obtain 
r d 4 k r 

J Tw ImI1 ^ ' ^ = ~ 2 sInn J ■ 



d A k r TTn n 2 /■ dV 3 ^ k\-k 2 



( 27T ) 4 ' 7 (2n)S2y/M* + {p - £) 2 2^ + p 2 

xe(pj,- |^l)e(|^- (23) 



with fc = yM 2 +p 2 — y M 2 + (p — k) 2 on the right-hand side. Comparison of expres- 
sions (|21|) and (|23|) shows that in general one cannot express (|2T|) in terms of (|23|) because 



of presence of the momentum-dependent factor l/JM 2 + p 2 in Eq.fl23|). However, if we 



approximate yM 2 + p 2 with a p-independent constant (say M*), which is an excellent ap- 
proximation for baryon densities not much larger than the saturation density, the following 
relationship emerges: 



L = -4M* / -0-I m Il(k o ,k). (24) 



Apart from the factor 1/y M 2 + p 2 whose momentum dependence in the present context 
does not play an important role, there is another momentum-dependent factor in expressions 



(^T|) and (|23|): 1/y M 2 + (p — k) 2 , whose presence is vital for correct low-density limit. In 
our derivation of expression (|22|) this factor appears automatically since we use the fully rela- 
tivistic nucleon propagators in the particle-hole loop. However, by using the non-relativistic 



limit (as, for example in Ref. |1|), this term would be replaced by 1/M*. This leads then to 
an incorrect low-density limit of the pion distribution, since the Sullivan expression @ for 
the free-nucleon is based on relativistic calculation and contains the momentum dependence 
of the mentioned term (with p = 0). 

Taking into account the more complicated diagrams means replacing II by Yl+UD®k 2 Il+ 

UDlk 2 YlDlk 2 Il H = 11(1 — _D°/c 2 n) _1 , where D® is the free-pion propagator. Also, when 

integrating over the pion momentum we have to take into account the constraint on the pion 
light-cone momentum fraction y = (ko + k^)/M. This can be imposed by an appropriate 
delta-function, which should be inserted in both expressions (0) and (|23|) . Performing the 
&3 integration and dividing by the total number of nucleons in unit volume we obtain the 
expression for the pion light-cone momentum distribution 

where (in the general case) the pion self-energy which contains both the delta-hole and 
particle- hole contributions is given by Eq. (|j). In II and D® the argument k is given by 



k = J kj_ + k 2 , where ks = My + uj. The integration region here is from zero to infinity for 
both variables, since conditions on incoming nucleon (below the Fermi sea) and outgoing 
nucleon (above the Fermi sea) are automatically satisfied by the correct calculation of the 
imaginary part of II, if the nucleons are treated as a mean field. 

Using the expression for the nucleon density p = 2p 3 F /3n 2 and introducing the longitu- 
dinal response-function through 
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, , N 1 4MM, / U(u,k) \ , \ 

R L (u, k) = — 2 * Im 1 ' , 26 

n g z NN F 7rNN (k) 2 \l- D^k 2 Il(uj,k) J 

expression fl25|) can be written in the usual form 

/7rM( } = /■ djW WJfc)2 / du 2 R f' k \ (27) 

lbir z p J J (uj z — k — m^Y 

Note that the form factor and the coupling have been removed from the imaginary part of 
the (iterated) pion self-energy in the definition (|26|) of the response function. 

We confirmed by numerical computation that Eqs. ( |ll|) and (|25|) give practically identical 
results (thus validating the approximation used to obtain eq.(|4])), if the pion self-energy 
contains only the particle-hole contribution (assuring the final state with only nucleon and 
not delta, corresponding to diagrams summed in Eq. (|TTD). Inclusion of the self-consistent 
delta-hole contribution to the pion self-energy takes into account also the diagrams where 
a delta replaces the nucleon in the final state. This corresponds to taking into account the 



contribution of diagrams summed in eq. (15), but also some other higher-order diagrams, 



which appear because the delta is itself dressed by pion-nucleon loops (where the pion is 
also dressed). Numerically this shows up as a very small difference (at values y < 0.2) of the 
distribution fl27D and the sum of eqs.(p|) and (15). The smallnes of this difference indicates 
that it is reasonable to expect that contributions of other higher-order diagrams left out 
from the calculation (for example, the pion-nucleon loop term of the nucleon self-energy) 
will not affect the pion distribution significantly. 

IV. DRELL-YAN SCATTERING 

The cross-section of the Drell-Yan process a + b — ► lT can be expressed as 

^ = ^M^) V' [M,{X2) + ql(xi)q ' {X2)] • (28) 

where s is the center-of-mass energy squared and the summation of products of quark and 
antiquark distribution functions is over flavors. The factor K{x\,X2) takes into account 
higher order QCD corrections and is of the order 1.5. The values of X\,x-2 are extracted 
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from experiment via the invariant mass of the lepton pair. We are interested in the ratio of 
the cross-sections for proton-nucleus and proton-deuteron scattering: 

2 da pA i 'dx\dx 2 



RA/d = 



(29) 



A da pd / dx\dx 2 ' 

where A denotes both the nucleus and its nucleon number. 

Since our in-medium computation is performed for the case of isospin symmetric nuclear 
medium, we specialize to the ratio of differential cross-sections on a nucleus consisting of 
equal numbers of protons and neutrons, and on a deuteron (for which we consider medium 
effects negligible) 



RA/d 





q P / A (x 2 ) + q] /A (x 2 ) 


+ q p f (x 1 ) 


'q P / A (x 2 ) + q n / A (x 2 )]} 




q p f (x 2 ) + q](x 2 ) 


+ ffixi) 


q p f (x 2 ) + q]{x 2 ) 


} 



(30) 



To relate the quark distribution of the bound nucleon to that of the free nucleon we start 
from (0), where for brevity we write out only nucleon terms: 



p i \ rv p,bare 

qj{x) = Zq^ 



f kx) + 1 £ ^f N / N {y) [qf aTe {x/y) + 2qf are (x/y) 

1 f 1 *L r /N( y) \ q f( x / y ) + 2 qf (x/y) 
o Jx y L 



(31) 



with 



Z = 1 



dyf N/N (y) = l - / dyr' N {y) 



(32) 



where the last equality in (|32|) is a requirement for flavor-charge conservation. Similarly, the 
starting expression for the quark distribution in a nuclear proton is 



, , ry p,bare 

Qf{ x ) = z Aq f 



» + \ f A ~f N/A (y) \q P f Mre (x/y) + 2qf bare (x/y) 
3 Jx y 1 1 1 



1 f 1 d A r IA {y) [ q j\ x / y ) + 2qf(x/y) 
o Jx y L 



(33) 



where we specialized to the case of isospin-symmetric nuclear matter. Adding and subtract- 
ing the expression (|31|) to expression (|33| ) and repeating the same procedure for the neutron 
(since only the sum of these two terms is relevant in isospin-symmetric medium) we get 
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f f (x) + %{z) = q%x) + qf are (x) f 1 f N ' N {y)dy - f ^f N/N (y)q P f Mre (x/y) 

jo Jx y 

-qf ar \x) J Q A f N/A (y)dy + l A ^f N/A (y)qf bare (x/y) + ip^n) 



+ 2 r^dy 
3 Jx y 



r /A (y)-r /N (y) qf(x/y) + qf(x/y) + qf(x/y) . (34) 



To proceed without approximation a fit of the bare structure functions q p,bare (x) would 
be required, based on expression fl3T|) and the experimentally extracted q V t{x). This is 
clearly beyond the scope of present work and in the following we present arguments that the 
expression used in Refs. |$,18,8| for the pion contribution to the change of parton distribution: 



5q p f (x)= f A ^-5r(y)q](x/y), (35) 
Jx y 



(where Sf n (y) = f^^iy) — f n ^ N (y) and isospin factors are not shown) is a good approx- 
imation for antiquark distributions, the change of which is probed by the Drell-Yan pair 
production. 

We argue that for antiquark distributions at small a; it is a good approximation to 
neglect the difference of the second and third term on the right-side of Eq. (0), as well as 
the difference of the fifth and fourth term. The reason is that f N ^ N (y) and f N ^ A (y) are very 
small for y < 0.3 (since the pion distributions are negligible for y > 0.7), thus if x < 0.3 (the 
region where antiquark distributions are significant), the zero lower limit of integrals can be 
safely shifted to x. Taking into account that antiquark distributions at small x behave as 1/x 
(and that x/y is also small in the y region of most significant contribution), cancellation of 
the considered terms follows. A much simpler argument in favor of expression (^) is based 
on the assumption that the nucleon's entire antiquark sea can be attributed to its virtual 
meson cloud which, however, has not been confirmed ||13|| . Since valence-quark distributions 
increase with decreasing x slower than 1/x in the small x region, the above approximation 
may be less good in that case, but the effect on the Drell-Yan process of a relatively small 
change in quark distribution is not significant. The sum of the nuclear proton and neutron 
(anti) quark distribution of flavor / thus becomes 
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f f (x) + q](x) = <f f (x) + q](x) + ~ t ^ \f^ A {y) - ^(y)} [qf(x/y) + qf(x/y) + qf(x/y) 

o J x y L J L 

(36) 



Performing a convolution on the nucleon part to take into account the Fermi motion of the 
nucleons we obtain 

" A dz 



q P / A (x) + qf A (x) 



/JLmO) [<ff(x/z) + q](x/z 

2 r A dy 



~*~3 Jx y 



F ,A {y) - r /N (y)] WAx/y) + qUx/y) + q^{x/y)] . (37) 



We note that the above expression satisfies flavor-charge conservation by construction 
and the momentum conservation sum-rule (with only pions included) 



f A ydy\r /A (y)-F /N (y) 



zdzfj, 



I. 



(38) 



'0 jo 

turns out to be satisfied with values of M* and cq based on nuclear matter models with 
accuracy better than 1% (see next section). Heavier mesons in general give much smaller 
contributions than the pion |L2] and would not affect significantly the sum-rule. 

One can treat the contribution of the delta in the same way, an important difference 
being that even in the isospin symmetric case the structure function of the delta appears in 
the structure function of the in-medium nucleon. 

To be consistent with mean-field treatment of nucleons which was used for calculating 
the pion's self-energy, for the nucleon light-cone momentum distribution /^ ermi we use the 



Fermi gas model with mean-field corrections for mass and energy, E{p) = y M% + p 2 + Co 
The expression, which takes into account the flux factor, is the one obtained by Birse |T9 



f w (z) 

J FermiK I 



4e 3 



(z — r/) 2 0(e — \z — r]\) 



(39) 



where e = p F /M,rj = (J M 2 + p 2 F + c )/M and Q(x) is the step function. 



V. NUMERICAL RESULTS AND DISCUSSION 



First we examine the effect of Pauli blocking and Fermi motion on the pion light-cone 
momentum distribution, f nA (y) by comparing in Fig. 3 the free-nucleon case (full line) to the 
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in-medium one, but the latter calculated with the free-pion propagator. The pion-nucleon- 
nucleon coupling has the usual value g n NN = 13.5. 

For the actual form of the form factor various parametrizations have been used, e.g. 
F(t) = (^zf) n with n = 1 (monopole) or n=2 (dipole), and also exponential forms. 
The corresponding values of the cut-off, A, can be related, see Ref. |H| . Previous studies 
|5| Jl^JT3|j20|j2l|| of observables related to antiquark distributions obtained cut-offs in the range 
0.6 — 0.8 GeV, when translated to monopole form, which is much smaller than the one used 
in the Bonn potential. We performed computations with a monopole form using A^at = 0.7 
and 0.8 GeV, and also exponential form with 1 GeV. While the effect of changing the cut-off 
A w nn is significant for both the free nucleon case and in the medium, their difference shows 
little sensitivity to this change. This is in contrast to Ref. || where a decrease in the cut-off 
produced a decrease in the enhancement too. In the following for all presented results we 
used a monopole form-factor with A^nn = 0.8 GeV, both in vacuum and in the medium. 

For the values of M* and Cq we turn to nuclear matter models. While the simplest form of 



the Walecka model J^j gives at saturation density M* ^0.7 GeV, more elaborate treatments 
tend to increase the effective mass to 0.8 — 0.85 GeV at saturation density, leading to better 
agreement with observables p4||2"5| . Since we use for the Fermi momentum pp = 0.256 GeV 



(i.e. slightly below the saturation density), we take M* = 0.85 GeV. Adopting a smaller 
value for the nucleon's effective mass decreases the in-medium pion enhancement, as noticed 
in Ref. || . The value of the energy shift cq should be such that a reasonable value of the 



Fermi energy is obtained ||19|| , thus we take cq = 0.04 GeV, giving Ep — M = — 11 MeV. 
For these values of M* and Cq the momentum-conservation sum-rule (|38|) is satisfied with 
accuracy better than 1% (for values of the g' parameters in the region from to 0.4), still 
leaving a little space for the effect of other mesons. 

In Fig. 3 we see a reduction due to Pauli blocking and practically no broadening for the 
effective mass of the nucleon M* = 0.85 GeV (short-dashed line). We note, however, that 
the decrease of the effective mass of the nucleon contributes to the decrease and narrowing 
of the pion distribution (as pointed out in Ref. [g]), as can be seen also from Fig. 3, where 
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the results for M* = 0.939 GeV and 0.7 GeV are also shown. 

For the in-medium calculation we consider first the contribution from diagrams with 
a nucleon in the final state. The results for fjJ A (y), Eq. ([□]), for different values of g' 
parameters are shown and compared to the free-nucleon case in Fig. 4. For g' NN we used 
the values 0, 0.4 and 0.6, while the other two parameters are kept constant, g' NA = g' AA = 
0.3. Increasing the g' parameters in general leads to reduction in the enhancement of pion 
distribution, although for small values of g' there can be some increase for certain y values, 
as a consequence of corresponding enhancement in some kinematical regions [|J. For the 
in-medium calculation we keep the mass difference of the delta and the nucleon the same as 
in free space. 

We also consider the contribution from the delta in the final state, f A A (y)- In connection 
with it we make some remarks on the 7riVA coupling. In the past the cut-off for the 7riVA 
coupling was taken the same as for the ttNN (using SU(6)). More recently there is evidence 
that the former coupling is softer than the latter from several sources: 

(i) The Juelich group has analyzed data for diffractive scattering p + p — > A ++ + X in the 
one-pion exchange approximation, 

(ii) Koepf et al. ]TB|| pointed out that the empirical u — d asymmetry can be explained in a 



meson cloud model only if A^na is softer than K^nn by about 100 MeV. 

In the present approach || we have extracted A^a from a fit to irN scattering phase- 
shift in the delta channel, using three parameters: mass of the delta, value of the coupling 
and cut-off, and obtained the value for the exponential form-factor of A n ^ A = 0.38 GeV 
(for dipole form-factor the value of the cut-off was 0.51 GeV). This value gives an excellent 
fit for the phase-shift for laboratory pion momentum up to 500 MeV. The tcNA coupling 
used, g-n-NA, was 20 GeV -1 (this implies an on-shell value of 14 GeV^ 1 , in good agreement 
with expectation based on the delta's width), and the (bare) delta mass of 1.27 GeV. We 
remark that a similarly soft pion-nucleon-delta form-factor was extracted from pion-nucleon 
scattering in Ref. P^ ], where also a self-consistent dressing of the pion and the delta was 
performed. Since the 7rA^A cut-off is rather small, we expect this contribution to be small. 
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This is indeed the case as one can see from Fig. 5. The pion distribution is nonzero only in the 
region of small y (y < 0.2), and while for g' — (here the g' AA plays the most prominent role) 
there is some enhancement over the free-nucleon case (the full line), its value is very small. 
To study the effect of the cut-off, we also performed computations with A^na = 0.5 GeV 
(which gives a very poor fit to the pion-nucleon phase-shift) and confirmed the absence of 
medium enhancement in that case too. We note that using the spectral-function of the 
dressed delta in the medium presents numerical difficulties, since it is nonzero in the region 
where the real part of the pion propagator takes very large values. To obtain stable results 
we neglected the delta contribution if its spectral function was smaller than 0.2 GeV -1 (this 
affects the spectral-function sum-rule negligibly). 

The advantage of the expression (^) relating the response function to the pion light- 
cone-momentum distribution is that it incorporates both contributions with a nucleon as well 
as delta in the final state, if the pion self-energy contains the particle-hole and the delta-hole 
contributions. Numerical evaluation confirms that for y > 0.2 these results agree with those 
based on expression (]TT|), while for small y there is a small enhancement due to diagrams with 
the delta in the final state. The slight difference between results based on expression ( |2"7D 
and the sum of contributions ([TT]) and flT5] ) can be attributed to contributions of higher order 
diagrams present in the self-consistent delta self-energy. The results of free-space (which is 
approximated by a low-density calculation with p = 0.1p , where po corresponds to nucleon 
density of pp = 0.256 GeV) and in-medium computations, using the same parameters as for 
Fig. 4, are shown in Fig. 6. 

As to the origin of the absence of significant in-medium enhancement of the pion cloud 
we point to the results of Ref. || for the pion spectral function. The basic feature of that 
model is the self-consistent treatment of the pion and the delta, with both real and imaginary 
parts of the self-energies taken self-consistently into account, thus assuring causality for self- 
energies and propagators. As a result the pion spectral function has a main maximum 
whose shift to smaller energies is much less pronounced than in previous, less elaborate 
treatments, an exception being only Ref. ||22|| , whose results show similar behavior to those 
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of Ref. ||. Since this shift forms the origin of the pion-cloud enhancement and comes from 
the (negative) real part of the delta-hole contribution, we can conclude that this is where 
the relevant difference between our pion's self-energy and that of previous models' rests. 
The sensitivity of the pion enhancement to this term requires its careful treatment, which 
is evidently missing in Refs. |3,r6|. We emphasize once more that in Ref. [[J the real part 



of the self-energy is always calculated from its imaginary part using a relevant dispersion 
relation. 

In the calculation of the structure functions of nuclear nucleons we do not take into 
account the delta contribution. The reason is that we do not know the (bare) delta structure 
functions which are needed as are the nucleon structure functions. The delta structure 
functions are necessary even in the case of the isospin symmetric nuclear matter. However, 
the above results for the pion light-cone distribution show that there can be no noticeable 
medium enhancement. 

We now turn to a discussion of the DY ratio, Eq. fl3Cf). In applying Eq. (|3~7| ) we use the 
quark distributions in the nucleon and the pion from Refs. [|26|| , and |27| , respectively. The 



data of Ref. cover a rather limited region of £2 values around x<± = 0.2. That means that 
in the expression (1371) the convolution in the first term on the right-hand side is relatively 



unimportant and that in the integral in the second term rather large y values are sampled, 
typically y > 0.3 and larger. Therefore one probes rather large 3-momenta in the response 
function D(k , k), Eq. (|7J). From k 2 > k\ and k 3 = —k Q + my one has \k\ > \k \ + my , i.e. 
typically k > 400MeV/c. As a consequence Pauli blocking (which is relevant for small k ) 
is not very effective; also the delta-hole component at k Q ^ is relatively unimportant since it 
involves still higher momenta, k > my + |/co,aI- This is rather different from the situation 
for the longitudinal spin response function probed in e.g. the quasi-free (p,n) reaction. 

Results for the ratio ([!(]) of Drell-Yan cross sections are shown in Fig. 7, for different 
values of the Migdal g' parameters. A general feature of all these plots is the decreasing 
trend for larger values of £2 (^2 > 0.3), leading to values below unity. The enhancement 
for smaller values of £2 even with no Migdal correction is quite modest. Values of the g' NN 



around 0.3 lead to practically no enhancement. 

Compared to Ref. || we still obtain a net medium DY enhancement by 10% for g' < OA 
in the region 0.1 < 22 < 0.3. This difference can be ascribed to the the fact that in |J an 
ad hoc Z-factor (see discussion below Eq. (fj)) occurs on right-hand side of Eq. ( |37| ) which 
is not present in our approach. In addition in || the assumed density dependent in-medium 
reduction of the cut-off A in the ttNN form factor appears the main mechanism to reduce 
the pionic enhancement. However we feel there is no justification for such a medium vertex 
renormalization. 

For direct comparison with measured results of Ref. we calculated the ratio of nuclear 
and nucleon cross-section for given values of X2 and the condition x\ > £2+0.2, corresponding 
to the experimental cut-off (here we assume that the factor K(xi 1 X2) is constant). In 
accordance with the above discussion we verified numerically that the ratio in the given x 2 
region is mainly sensitive to the pion distribution in the region y > 0.3. Only the region of 
X2 around 0.3 shows some sensitivity to smaller y values. This means that the drop of the 
ratio below one for x 2 around 0.05 (if the effect is real), probably cannot be explained by a 
change of the pion distribution. 

Based on our results we can conclude that the treatment of the delta, although seemingly 
of secondary importance because of the kinematical region of the delta-hole contribution, 
does play a significant role. Including both the real and imaginary parts of its self-energy in 
a self-consistent way eliminates a large part of the pion-distribution in-medium enhancement 
and leads to agreement with experimental results on Drell-Yan scattering for modest values 
of the g' parameters. The reduction in the in-medium pion enhancement is achieved without 
intruducing a medium reduction of the pion-nucleon-nucleon vertex cut-off parameter or a 
renormalization factor to suppress the pion-cloud contribution. 

Finally we note that in the present paper we have restricted ourselves to isospin symmet- 
ric nuclear matter. On the other hand some of the data involve nuclei like Fe and W, which 
have a neutron excess. It is an interesting question whether there are additional medium 
effects to the asymmetry u(x) — d(x) in case of neutron excess. We plan to address this 
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question in the future. 
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FIGURES 

FIG. 1. Types of diagrams involving the pion self-energy and dressed propagator (dashed line), 
used for calculating the pion distribution in the medium. The double line denotes either a nucleon 
or delta. 

FIG. 2. The particle-hole pion self-energy (a), and the pion emission part of the diagram 
showing the deep- inelastic scattering off the nucleon's pion cloud (b). 

FIG. 3. Effect of Pauli blocking and Fermi motion on the pion distribution. The full line 
corresponds to to the free nucleon, the short-dashed line to nuclear matter with pp = 0.256 GeV 
and M* = 0.85 GeV, but with free-pion propagator. The long-dashed line corresponds to M* = 0.7 
GeV, while the dot-dashed line is for M* = 0.939 GeV. 

FIG. 4. Pion distribution in the nuclear medium with Fermi momentum pp = 0.256 GeV 
and for the free nucleon. Full line is for the free nucleon, long dashed line for medium with 
9nn = 0' 9na = 9aa = short dashed for g' NN = 0.4, g' NA = g' AA = 0.3 and dot-dashed for 
9nn = °- 6 > 9na = 9aa = °- 3 - 

FIG. 5. Delta contribution to the pion distribution in the nuclear medium with Fermi momen- 
tum pf = 0.256 GeV and for the free nucleon. Full line is for the free nucleon, long dashed line for 
medium with g' NN = g' NA = g' AA = 0, short dashed for g' NN = g' NA = g' AA = 0.3 and dot-dashed 
for g' NN = g' NA = g' AA = 0.5. The effective mass of the delta in the medium is M A = 1.18 GeV 
(assuring the same effective mass difference for the delta and nucleon in the medium as in the free 
space), and the 7T./VA vertex cut-off A^na = 0.38 GeV. 

FIG. 6. Pion distribution in the nuclear medium, based on computation using the re- 



sponse-function approach, expression (27). Full line is for the low-density case with 
PF = 0.119 GeV, corresponding to 1/10 of the density used for in-medium calculations. 
Long dashed line is for medium with g' NN = g' NA = g' AA = 0, short dashed for 
9nn = oa >9na = 9aa = °- 3 and dot-dashed for g' NN = 0.6,g' NA = g' AA = 0.3. 
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FIG. 7. Ratio ( j30|) of the Drell-Yan cross-sections. Full line is for x\ = 0.3, dashed line for 
x\ = 0.4 and dot-dashed line for x<i = 0.5. The Migdal parameters are: a) g' NN = g' NA = g' AA = 0, 

b ) 9nn = oa ,9na = 9aa = °- 3 - 

FIG. 8. Experimental results from Ref. Q compared to our calculation for different g' values. 
Full line is for g' NN = g' NA = g' AA = 0.4, dashed line for g' NN = 0.4, g' NA = g' AA = 0.3, dot-dashed 
line for g' NN = 0.6, g' NA = g' AA = 0.3. 
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Fig. 7b 
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